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Abstract. Quantitative analyses of low-mass hydrogen-deficient (super-)giant stars - so-called extreme helium stars - to date 
face two major difficulties. First, theory fails to reproduce the observed helium lines in their entirety, wings and line cores. 
Second, a general mismatch exists for effective temperatures derived from ionization equilibria and from spectral energy distri- 
butions. Here, we demonstrate how the issue can be resolved using state-of-the-art non-LTE line-formation for these chemically 
peculiar objects. Two unique high-gravity B-type objects are discussed in detail, the pulsating variable V652 Her and the metal- 
poor star HD 144941. In the first case atmospheric parameters from published LTE analyses are largely recovered, in the other 
a systematic offset is found. Hydrogen abundances are systematically smaller than previously reported, by up to a factor ~2. 
Extreme helium stars turn out to be important testbeds for non-LTE model atoms for helium. Improved non-LTE computations 
show that analyses assuming LTE or based on older non-LTE model atoms can predict equivalent widths, for the He 1 10 830 A 
transition in particular, in error by up to a factor ~3. 
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^ 1. Introduction 

Extreme helium stars (EHes) are a rare class of hydrogen- 
deficient objects showing spectral characteristics of early A 
to late O-type (super-)giants. Their chemical composition is 
dominated by fusion products from the CNO-cycle and/or 
the 3a-process. Therefore, EHes provide important clues for 
the study of nuclear astrophysics. Most of the two dozen 
known EHes could be explained by post-AGB evolution, link- 
ing R Cr B stars to Wolf-Rayet type central stars of planetary 
nebulae, see Heber J1986> and JefFery J 19961 1 for reviews. 

V652Her and HD 144941 are unique among the class 
members in several aspects. The carbon, nitrogen and oxy- 
gen abundances of both stars are consistent with being pro- 
cessed through the CNO cycle, while all other EHes display 
C-rich material at the surface, indicating that Sff-processing 
has occurred in addition. Both stars have gravities too large 
for post-AGB evolution. Saio & Jeffery (200(fi suggest that 
V652 Her may be the result of a He-i-He white dwarf binary 
merger. HD 144941 is also outstanding because of its strong 
metal-deficiency, larger than observed in any other EHe star. 

Because of their peculiar chemical composition EHes are 
important testbeds for stellar atmosphere modelling. Since he- 
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lium is by far the most abundant element, the He line spectrum 
can be studied in more detail than in any other type of star, e.g. 
in HD 144941 all forbidden transitions of He i can be measured. 
Hydrogen is deficient by a factor of 100 or more and therefore 
the Balmer lines are very weak, unlike in any other B-type star. 

Several LTE spectral analyses encountered two difficulties 
to be most troublesome: i) synthetic spectra have so far not suc- 
ceeded in matching the observed helium lines in their entirety, 
and ii) spectroscopic and spectrophotometric temperatures dif- 
fer systematically. As inadequacies in the basic parameter de- 
termination can potentially hamper any further interpretation, 
the issue needs to be resolved. In the following we investigate 
the steps that need to be taken to improve the modelling, and 
the consequences of these for our understanding of EHe stars. 

2. Model calculations 

Our model calculations are carried out in a hybrid approach, 
as commonly used for analyses of analogous B-type stars 
with normal chemical composition. First, hydrostatic, plane- 
parallel and line-blanketed (via an opacity sampling technique) 
LTE model atmospheres are computed with the Atlas 12 code 
(Kurucz 1996). Note that we have replaced the photoioniza- 
tion data for He i levels with principal quantum number n = 2 
as used by Kurucz with data from the Opacity Project (Fernley 
et al. .l987j . In particular the cross-sections for the 2p level 
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Fig. 1. Comparison of the normalised spectra of V652 Her (lower) and the metal-poor hydrogen-deficient star HD 144941 (upper 
histogram) with our non-LTE spectrum synthesis for helium and hydrogen (full red line): excellent agreement for the entire line 
profiles - wings and line cores - is found. The He and H features have been labelled. Forbidden components of He i lines missing 
in our modelling are indicated by short vertical marks. The red wing of He 1 4922 A in HD 144941 may be affected by an artifact 
in the data reduction. 



are increased by a factor ~2 at threshold, thus improving the 
fits of computed energy distributions with observation. Stellar 
parameters and elemental abundances from Jeffery et al. (.2001. 
JWPOl) and Harrison & Jeffery ( fT997l HJ97) are adopted for 
the initial models of the stars. Details on the observations can 
also be found there, which are complemented by further data 
from Jeffery et al. ( IT^^ . 

Then, non-LTE line formation is performed on the resulting 
model stratifications. The coupled radiative transfer and statis- 
tical equilibrium equations are solved and spectrum synthesis 
with refined line-broadening theories is performed using Detail 

Both 



Table 1. Photospheric parameters of the sample stars 



and Surface (Giddings 1981 Butler & Giddings [19851 
codes have undergone major revisions and improvements over 
the past few years. State-of-the-art model atoms for hydrogen 
(Przybilla & Buder 2004 PB04) and heUum (Przvbilla 2005. 
P05) are utilised, which differ from older model atoms mostly 
by use of improved collision data from ab-initio computations 
for electron impact excitation processes. For He i these are the 
data from Bray et al. (2000) and Sawey & Berrington ( 1993l, 
and for Hi the data from computations by Butler, as sum- 
marised in PB04 - see this work and P05 for further de- 
tails on the modelling. For comparison, additional calcula- 
tions are made using the helium model atom of Husfeld et al. 
i tT^ HBHD89). The theory of Barnard et al. C1969 ) which is 
discussed further by Auer & Mihalas | |1972> is utilised for a 
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22000±500 
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4+1 
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realistic description of line broadening, supplemented by data 



from Dimitrijevic & Sahal-Brechot (I1990> . 

The spectrum synthesis is compared to observation in order 
to derive improved values for the stellar parameters and hydro- 
gen and helium abundances, giving 'best fits' in an iterative 
approach. Effective temperatures Teff are determined from the 
He i/ii non-LTE ionization equilibrium and the Stark-broadened 
He I lines act as surface gravity indicators. Stellar parameters of 
the final models (with estimated uncertainties) are summarised 
in Table including microturbulence ^. For V652 Her the at- 
mospheric parameters agree very well with those found by 
JWPOl, except for the hydrogen abundance nu, which is re- 
duced by a factor ~2 (see Sect. 13. 2> . For HD 144941, however, 
the resulting atmospheric parameters differ significantly from 
previous work of HJ97, implying a reduction in T^ff by 1 200 K 
and an increase in surface gravity by a factor ~2. The reduction 
in the hydrogen abundance is less pronounced in this star. 
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Fig. 2. Departure coefficients bi for Hei in the model of 
V652 Her as a function of Rosseland optical depth Tjoss for the 
n-l and 2 and selected n-3 and 4 levels from computations 
using the P05 and the HBHD89 model atoms (upper panel, 
full and dotted lines, respectively). The Hen ground state is 
also indicated (dashed line). In the lower panel ratios of line 
source function Si to Planck function By for several Hei lines 
are displayed. Line formation depths are indicated. The con- 
siderable differences in the populations of the 2s state in the 
two model atoms translate to a discrepant S i/By ratio for the 
Hei 1.083 yum transition, giving a substantially weakened line 
for the modern model atom, while computations based on the 
old model deviate only little from detailed balance. 



3. Testbed for non-LTE spectrum synthesis 

3.1. Optical helium line spectrum 

As can be seen from Fig. \l\ excellent agreement of the non- 
LTE spectrum synthesis with the observed line profiles of the 
He features - wings and cores (and forbidden components) - is 
finally obtained in our approach. Similar agreement of theory 
with observation in the visual spectral range is obtained using 
the HBHD89 model atom. The great improvement achieved be- 
comes obvious when comparing Fig. [2 to Fig. 2 of HJ97 and 
Fig. 5 of JWPOl. This resolves one of the most persistent in- 
consistencies in quantitative analyses of EHes. Only a few pre- 
dicted forbidden components (Beauchamp & Wesemael (199811 
are missing because the broadening data are unavailable to us. 

The He i lines in the visual experience non-LTE strength- 
ening, facilitated by the non-LTE overpopulation of the n-2 
states relative to the levels of higher principal quantum number, 
see Fig. |2l There, the runs of non-LTE departure coefficients 
bi - tiiln* (the n, and n* being the non-LTE and LTE popula- 
tions of level /, respectively) and line source functions Si (rel- 
ative to the Planck function By) are displayed. The non-LTE 
overpopulation occurs because of recombinations to levels of 
He I at higher excitation energies and subsequent de-excitation 
via downward cascades to the (pseudo-)metastable 2s states 
(the singlet resonance lines are close to detailed balance). The 
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Fig. 3. Same as Fig. |2] but for Hi. The are labelled by 
the principal quantum number. Non-LTE (full line) and LTE 
profiles (dashed line) of Hy6 are compared in the inset. Line 
source-functions are displayed for Balmer lines and Paschen 
and Brackett lines in the J- and /T-band. 

singlet lines are in general subject to larger non-LTE effects 
than the triplet lines. Note that the level populations deviate by 
only a few percent from detailed equilibrium at the formation 
depths of the continua, showing that the assumption of LTE for 
the model atmosphere computations is appropriate. 

3.2. Optical hydrogen line spectrum 

In analogy to helium hydrogen can also be supposed to expe- 
rience non-LTE effects. Departure coefficients for several en- 
ergy levels and line source-functions for selected H i lines are 
displayed in Fig. |3l The Hi ground state is strongly depop- 
ulated, however at depths irrelevant for the formation of the 
Lyman continuum. Consequently, non-LTE effects on H are 
also unimportant for the atmospheric structure calculations. A 
small overpopulation of the first excited level occurs at line- 
formation depths. All other levels of higher n are closely cou- 
pled to the continuum (which is in detailed balance) at line- 
formation depths. This results in a non-LTE strengthening of 
the Balmer lines, with non-LTE and LTE equivalent widths for 
e.g. Wfi differing by more than 30%. A reduction of H abun- 
dances by a factor up to ~2 is indicated relative to previous 
LTE studies (for V652 Her the lines are formed on the flat part 
of the curve of growth). 

3.3. Infrared line spectra 

Spectral lines in the Rayleigh-Jeans tail of the spectral energy 
distribution, e.g. in the near-IR for EHes, can experience ampli- 
fied non-LTE effects (see e.g. PB04). Analyses in the near-IR 
range are highly useful for constraining the atomic data input 
for the non-LTE computations. He i in the present case. 

Non-LTE weakening is indicated from Fig.|3for the 2s^S- 
2p transition in the 7-band, and an emission feature in the 
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Fig. 4. Comparison of the 7-band spectrum of V652 Her (his- 
togram) with three model predictions: LTE (dashed line) 
and two non-LTE computations using the model atom of 
HBHD89 (dotted line) and the improved model atom of 
P05 (full red line). In contrast to the other transitions, He: 
/ll.083yum experiences strong non-LTE weakening. An un- 
broadened synthetic spectrum is displayed at the bottom in or- 
der to illustrate what can be expected from high-resolution ob- 
servation. 



case of the 2s ^S-2p 'f ° transition in the AT-band. Here, the 
upper levels are overpopulated relative to the lower states. The 
He I 1.083/im transition turns out to be highly sensitive to the 
details of the model calculations, in particular to the use of ac- 
curate photoionization cross-sections for the 2s level, and 
the account of line blocking. As is the case of early-type main 
sequence stars (P05), large differences to computations based 
on older He i model atoms such as those of HBHD89 are found. 
The state-of-the-art model atom predicts the equivalent width 
Wa of He I /iL083;um to be a factor ~3 smaller than the old 
model atom, which indicates only small deviations from de- 
tailed balance. 

The significance of the predictions based on the new model 
atom is demonstrated in Fig. 0] where a comparison is made 
with the only published spectrum for an EHe star in the 
near-IR. Good agreement of the model computations with the 
intermediate-resolution /-band observation (JWPOl) is found, 
containing in addition many strong He i lines originating from 
n-3 states. Note however, that the line broadening is dom- 
inated by instrumental effects. High-resolution observations 
would be highly desirable in order to constrain the modelling 
even further. The analogous 2 s ^S-2p 'P° transition in the K- 
band at 2.059;um is predicted to be entirely in emission, with 
Wx- - 200 mA for V652 Her The non-LTE effects are smaller 
for HD 144941 because of higher atmospheric densities, imply- 
ing a lower ionization of the atmospheric heUum. Finally, we 
find that the hydrogen Paschen and Brackett series are much 
closer to detailed equilibrium than the Balmer lines. 

4. Spectral energy distribution 

The second major difficulty in the modelling of EHes is a gen- 
eral mismatch of effective temperatures derived from ionization 
equilibria and spectrophotometry. We compare the Atlas 12 
model fluxes for stellar parameters derived from the He i/n ion- 




log A (A) 

Fig. 5. Fits of theoretical energy distributions (full red lines) to 
lUE, visual and IR spectrophotometry for V652 Her (near max- 
imum radius, lower jagged curve and dots) and HD 144941 (up- 
per curve). The observed fluxes are dereddenedby the amounts 
indicated and the HD 144941 data are shifted vertically by -1-0.3 
units. The models are normalised to the observed Johnson V 
magnitudes. Locations for the bound-free edges of the n-2 
levels of He i are indicated. 

ization equilibria (see Tabled with lUE spectrophotometry 
and visual and near-IR photometry in Fig.|5] Excellent agree- 
ment is found for V652Her, and reasonable agreement for 
HD 144941 when interstellar reddening is accounted for 

The fit for HD 144941 may be improved by allowing for 
changes of the stellar metallicity and/or the effective tempera- 
ture. Also, a non-standard reddening law cannot be ruled out. 
A completely satisfactory solution for this metal-poor object 
is beyond the scope of the present paper, as metal abundances 
have to be determined, and T^^ should be verified from non- 
LTE ionization equilibria of metals. However, model atoms for 
several of the strategic species are unavailable at present. 

5. Conclusions 

The present work demonstrates how a hybrid non-LTE ap- 
proach using line-blanketed LTE model atmospheres and non- 
LTE spectrum synthesis succeeds in solving the most persis- 
tent problems in the quantitative spectroscopy of EHes. Test 
computations for two objects imply that stellar parameters from 
classical analyses of EHes may be subject to systematic modifi- 
cation in the hybrid non-LTE approach. Lower hydrogen abun- 
dances than reported previously are found, by up to a factor ~2. 

The atmospheres of the chemically peculiar EHes allow 
non-LTE model atoms in an environment complementary to 
normal stars to be tested. The present work verifies the inade- 
quacy of LTE and older non-LTE model atoms for helium. Only 
a recently improved model atom (P05) reproduces observation 
for all available spectroscopic indicators, including the Hei 
/II. 083 yum feature in particular, which changes by a factor ~3 
in equivalent width while the lines in the visual remain practi- 
cally unaffected. High-resolution spectroscopy in the near-UV 
and near-IR will be required to put even tighter constraints on 
the atomic model for one of the most basic elements. Finally, 
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an improved understanding of EHes in the context of nuclear 
astrophysics will require non-LTE studies of metal abundances. 
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